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£2 ■ Abstract 

We analyze the rare kaon decays Kjj — > ir°i>9, K + — > ir + i/i>, Kl —* 7r°e + e _ and i?£ — > 
/i + /U _ in conjunction with the CP violating ratio e'/e in a general class of supersymmetric 
models in which Z- and magnetic-penguin contributions can be substantially larger than 
in the Standard Model. We point out that radiative effects relate the double left-right 
mass insertion to the single left-left one, and that the phenomenological constraints on the 
latter reflect into a stringent bound on the supersymmetric contribution to the Z penguin. 
Using this bound, and those coming from recent data on e'/e, we find BK(Kl — > tt^uv) % 
1.2-10 -10 , BR(K+ -> ir+vv) < 1.7- 10 -10 , BR(K L vr e + e-) d ir Z 2.0-10 -11 , assuming the 
usual determination of the CKM parameters and neglecting the possibility of cancellations 
among different supersymmetric effects in e'/e. Larger values are possible, in principle, 
but rather unlikely. We stress the importance of a measurement of these three branching 
ratios, together with improved data and improved theory of e'/e, in order to shed light 
on the realization of various supersymmetric scenarios. We reemphasize that the most 
natural enhancement of e'/e, within supersymmetric models, comes from chromomagnetic 
penguins and show that in this case sizable enhancements of BR{Kl — > 7r°e + e~)di r can 
also be expected. 



1 Introduction 



Flavour-Changing Neutral Current (FCNC) processes provide a powerful tool for 
testing the Standard Model and the physics beyond it. Of particular interest are the 
rare kaon decays Kl — > tt°uu, K + — ► n + uu and — > 7r°e + e~ which are governed 
by Z-penguin diagrams. The latter diagrams play also a substantial role in the CP 
violating ratio e'/e. The most recent experimental results for this ratio, 

Re( £ 7 £ )=( (28 - 0±4 - 1) - 10 " 4 (KTeV) W (1) 
{ (18.5 ± 7.3) • 10~ 4 (NA48) 

are in the ball park of the earlier result of the NA31 collaboration at CERN, (23.0 ± 
6.5) • 10 -4 [[J, and substantially higher than the value of E731 at Fermilab, (7.4 ± 
5.9) • 10~ 4 [|J. The grand average (according to the PDG recipe) including NA31, 
E731, KTeV and NA48 results, reads 

Re(e'/e) = (21.2 ± 4.6) ■ 10 -4 , (2) 

very close to the NA31 result but with a smaller error. The error should be further 
reduced once complete data from both collaborations will be analyzed. It is also 
of great interest to see what value for e'/e will be measured by KLOE at Frascati, 
which uses a different experimental technique than KTeV and NA48. 

The estimates of e'/e within the Standard Model (SM) are generally below the 
data but in view of large theoretical uncertainties stemming from hadronic matrix 
elements one cannot firmly conclude that the data on e'/e imply new physics [^. 
[7], H . On the other hand the apparent discrepancy between the SM estimates and 
the data invites for speculations about non-standard contributions to e'/e. Indeed 
the KTeV result prompted several recent analyses of e'/e within various extensions 



of the Standard Model (see e.g. [[M]]) and particularly within supersymmetry [11, 12 
Unfortunately these extensions have many parameters and if only e'/e is considered 
the analyses are not very conclusive. 

The approach we want to pursue in the present paper is different: we will adopt 
a model-independent point of view within a generic supersymmetric extension of the 
Standard Model with minimal particle content, and study what are the implications 
of a supersymmetric e'/ e for the rare decays. To do so we will use the mass-insertion 
approximation [T3. Despite the presence of a large number of parameters within 
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this framework, only a few of them are allowed to contribute substantially to e'/e. 
Phenomenological constraints, coming mainly from AS = 2 transitions |14j, make 
the contribution of most of them to AS = 1 amplitudes very small compared to 
the Standard Model one. The only parameters which survive are the left-right mass 
insertions contributing to the Wilson coefficients of Z- and magnetic-penguin oper- 
ators. As we will discuss below, the reason for this simplification is a dimensional 
one: these are the only two classes of operators of dimension less than six contribut- 
ing to e'/e. Supposing that the enhancement of the Wilson coefficients of either of 
these two (or both) type of operators is responsible for the observed value of e'/e, a 
corresponding effect in the rare decays should be observed. In what follows we will 
analyze in detail the relations between the size of the effect in e'/ e and those in the 
rare decays. 

The same kind of logic was already followed by two of us in JT5fl. There, this kind 



of analysis was carried through under the assumption that the dominant effect in 
AS* = 1 transitions was only an enhanced sdZ vertex. This analysis was motivated 



by an observation of another two of us [16] that the branching ratios of rare kaon 



decays could be considerably enhanced, in a generic supersymmetric model, by large 
contributions to the effective sdZ vertex due to a double left-right mass insertion. 
This double mass insertion had not been included in earlier analyses of rare kaon 



decays in supersymmetry [17], [18[ . In the latter papers only single mass insertions 
were taken into account, leading to modest enhancements of rare-decay branching 
ratios, up to factors 2-3 at most, as opposed to the possible enhancement of more than 
one order of magnitude allowed by the double mass insertion [16|]. The conclusion 



of the analysis in [15] was that the data on e'/e may constrain considerably the 
double left-right mass insertion and the corresponding enhancement of the rare- 
decays branching ratios. 

In the present paper we will improve the analysis in (jHJ with the aim to answer 
the following questions: 



Can the large double mass insertions suggested in |T6[ be further constrained? 
As we will see this is indeed the case. 



What is the impact of these new constraints on the analysis in |T5 ] 



? 



What is the impact on this analysis of contributions from chromomagnetic and 
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7-magnetic penguins to e'/e and K L — > n°e + e respectively? 

As we mentioned above, in generic supersymmetric theories a sizable contribution 
to e'/e could also be generated by the chromomagnetic-dipole operator. Actually, 
within supersymmetric models with approximate flavor symmetries, the latter mech- 
anism seems to be more natural than a strong enhancement of the sdZ vertex 



TjJ . Interestingly, if the Wilson coefficient of the chromomagnetic-dipole opera- 
tor gets enhanced, one should also expect a sizable effect in the branching ratio of 
Kl — > 7r°e + e~, due to the 7-magnetic penguin. In fact their Wilson coefficients 
receive contributions from the same type of mass insertion. 

The paper is organized as follows: In Section 2 we identify the dominant SUSY 
contributions to | AS\ = 1 amplitudes as those of dimension less than six. In Section 3 
we summarize the effective Hamiltonian for | AS\ = 1 transitions concentrating on the 
operators of dimension four (effective sdZ vertex) and five (magnetic penguins) and 
their corresponding Wilson coefficients. Here we introduce three effective couplings 
which characterize the supersymmetric contributions to the Wilson coefficients of 
these operators: A t for the Z penguin and for the magnetic ones. In Section 
4 we collect the basic formulae for e'/e and rare kaon decays in terms of these 
effective couplings. In particular we calculate the magnetic contributions to e' / e and 
K L — > ir°e + e~. In Section 5 we analyze indirect bounds on the effective couplings. 
The main result of this section is an improved upper bound on |A t | coming from 
renormalization group considerations. In Section 6 we present a detailed numerical 
analysis of rare kaon decays taking into account the recent data on e'/ e, the present 
information on the short distance contribution to BR(f^ — > /i + /U~) and the bounds 
on effective couplings derived in Section 5. Analyzing various scenarios we calculate 
upper limits on BR{K L 7c°e + e-) dh , BR{K L -> vrVz/) and BR(K + ir+up). We 
present a summary and our conclusions in Section 7. 



2 SUSY contributions to |A5| = 1 amplitudes 

In the Standard Model FCNC amplitudes are generated only at the quantum level. 
The same remains true also in low-energy supersymmetric models with unbroken R 
parity, minimal particle content and generic flavour couplings. The flavour structure 
of a generic SUSY model is quite complicated and a convenient way to parametrize 
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the various flavour-mixing terms is provided by the so-called mass-insertion approx- 



imation |fL3|j . This consists in choosing a simple basis for the gauge interactions 
and, in that basis, to perform a perturbative expansion of the squark mass matri- 
ces around their diagonal. In the following we will employ a squark basis where all 
quark-squark-gaugino vertices involving down-type quarks are flavor diagonal. 

In the case of |Aif>| = 1 transitions we can distinguish between two large classes 
of one-loop diagrams: 

• Box diagrams. These are present both in | AS\ = 1 and \AS\ = 2 amplitudes. In 
both cases the integration of the heavy degrees of freedom, associated with the 
superpartners, necessarily leads to effective four-quark operators of dimension 
six. The Wilson coefficients of these operators are therefore suppressed by two 
powers of a supersymmetry-breaking scale, that we generically denote by Ms- 
Here 1/M| plays a role similar to 1/Myy in the SM case. 

Since any mass- insert ion carries at most \AS\ = 1, the leading contribution 
to | AS | = 2 transitions starts at second order in this expansion. Denoting 
by 5 the generic ratio of off-diagonal terms over diagonal ones in the squark 
mass matrices, the coupling of \AS\ = 2 effective operators turns out to be of 
0(5 2 /M|). This has to be compared with the dominant SM coupling that is 
of 0(\f/Myy), where X t = V t *V t d- If we then impose that the supersymmetric 
contribution to \AS\ = 2 amplitudes is at most of the order of the SM one, we 
find 

6 /M s < X t /M w . (3) 

In the case of \AS\ = 1 amplitudes, the leading supersymmetric contribution 
starts already at first order in 5, similarly to the SM one that is linear in 
Xt. However, the dimensional suppression factor is always 1/Mj in the SUSY 
case and in the SM one. Therefore, if Ms 3> Mw, the constraint (§) 

implies that the supersymmetric contribution to \AS\ = 1 box diagrams is 
suppressed with respect to the SM one. This naive argument is confirmed by 
the detailed analysis of |TJ]], where it has been shown that |A5| = 2 constraints 



always dominate over \AS\ = 1 ones, as long as we consider only dimension-six 
operators generated by box diagrams with gluino exchange. 

Penguin diagrams. At the one-loop level this kind of diagrams is present only in 
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\AS\ = 1 amplitudes. Effective operators with lowest dimension generated by 
photon and gluon penguins are the so-called "magnetic" operators of dimension 
five. The coupling of these operators is of 0(8 /Ms) and therefore potentially 
competing with the SM contributions even if we impose the bound (Q). This 
naive conclusion is again confirmed by detailed analyses of gluino mediated am- 



plitudes [[L4j]. In this context it is found that only the chromomagnetic operator, 
induced by c?i,(fl) — SR(L) mixing, could lead to sizable (> 10~ 3 ) contributions 
to e'/e without violating any constraints from e. 

A different situation occurs in the case of Z-penguin diagrams, where the break- 
ing of SU (2) l allows to build an effective dimension-four operator of the type 
SLl^dLZ^. Denoting by Cz the dimensionless coupling of this operator, the 
integration of the heavy Z field leads to an effective four-fermion operator pro- 
portional to Cz/M\ without any explicit I/M5 suppression. This potential 
enhancement is partially compensated by the fact that the leading contribu- 
tion to Cz arises only at second order in the mass-insertion ||16|| . However, the 
absence of any I/M5 suppression makes this term particularly interesting both 
for rare decays [HJ and e' / e [f3J . 

Given the above considerations, in the following we will restrict our attention only to 
the dominant SUSY effects in |A5| = 1 amplitudes: those generated by the "mag- 
netic" dimension-five operators, induced by gluino exchange, and those generated by 
the sdZ vertex mediated by chargino exchange. Interestingly, under this assumption 
only the off-diagonal left-right entries of squark mass matrices are involved, in par- 
ticular the dL(R) — $r(l) mixing for the magnetic operators and the u^'^ — Ir one for 
the sdZ vertex. 

What we will not consider are the gluino and the chargino contributions to irre- 



ducible dimension-six operators. The former have been explicitly calculated in ||14| 
and found to be negligible, the latter are suppressed by (9(M^/M|) with respect to 
the corresponding contributions mediated by the sdZ vertex. To control the accu- 
racy of our approximation, we have explicitly checked that the impact of these terms 
is below 10%, with respect to the dominant ones, for squark/gaugino masses above 
~ 300 GeV. Finally, we will completely ignore the neutralino contributions which 
are known to be negligible due to the smallness of both electroweak and down-type 



Yukawa couplings [18 
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Since a large sdZ vertex is already present in the SM, the corresponding SUSY 
corrections can be easily incorporated without modifying the structure of the SM 
| AS 1 1 = 1 effective Hamiltonian. On the other hand, the dimension-five operators, 
neglected within the SM, require an adequate treatment and will be discussed in 
detail below. 



3 Effective Hamiltonian 

3.1 Operators and Wilson Coefficients 

On the basis of the discussion in the previous section, we introduce here the effective 
Hamiltonian containing all the relevant operators of dimension smaller than six. The 
only dimension-four operator of interest is the one given by the sdZ vertex: 

= -^ M ^ z ^ z ^ + h ' c ' • (4) 

where 

Z ds = \ t C (x t ) + X t H (x qx ) . (5) 

Here the first term on the r.h.s is the Standard Model contribution (evaluated in the 't 
Hooft-Feynman gauge) and the second one represents the dominant supersymmetric 
effect. The couplings A t and X t are defined by 

A* = V*V td , A t = (6 u LR ) 23 (6 u LR y i3 , (6) 

where Vij are the elements in the CKM matrix and, denoting by Mp D ^ the squark 
mass matrices, 

Explicit expressions for the functions C and H will be given below. 

The magnetic operators of dimension five appear in the effective Hamiltonian in 
the following way: 

Hit = + C-Q- + C 9 + Q+ + C-Q-) + h.c. , (8) 

where we have chosen the following operator basis: 

Q± = ^(s L a^F, u d R ±s R a^F^d L ) , (9) 

Qt = T l~ 2 (s L ^t a G%d R ±s R a^t a G; u d L ) , (10) 
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Full expressions for the Wilson coefficients generated by gluino exchange at the 



SUSY scale can be found in JTJJ. We are interested here only in the contributions 
proportional to l/m g , which are given by 

Tia s {mg) r. 



m g 



I^) n ±(^)JfoW 



(11) 
(12) 



where the are defined in (|7|) and the functions F and G are given in (pD| ) and 
In the (Q^ ,Q~) basis, the leading order anomalous dimension matrix reads 



/ 8/3 \ 
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(13) 



\ 32/3 4/3 J 

Therefore, integrating out SUSY particles at the scale m g > m t , one has 



Cf(m c ) = V C±(mg), 



C±(m § ) + 8 (1 - tT 1 ) Cf{m 9 ) 



where 



V 



<*s(mg)\ 21 ( a s {m t )\ 23 j a s (m b ) 



2_ 

2.1 



(14) 
(15) 

\a s {m t ) ) \a s (m b ) J \a s (m c ) J 
The dimension- five operators in (§) in principle mix also with Q 2 , the leading 
dimension-six operator of the SM \AS\ = 1 effective Hamiltonian (see e.g. |T9| ). 
However, the effect of this mixing can be neglected as long as we are interested in 
large enhancements of the Wilson coefficients of the dimension-five operators with 
respect to the SM case (more than one order of magnitude in the imaginary parts, as 
suggested in ||1T[). Therefore, as first approximation, in the following we will neglect 



the mixing of with Q 2 . 



3.2 Basic Functions 

The basic functions relevant for our analysis are 



Bn(x) 



X 



+ 



xln(x) 
1 — x ' (x — l) 2 
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(17) 



C (x) 



Hn(x) 



x 



X[X 



X 



6 , 3x + 2 



x — 1 (x — l) 2 
3 — 6x 2 + 3x + 2 + 6xln(x)) 



48(1 -xY 



Fo{x) 
G (x) 



4x(l + 4x-5x 2 + 4x ln(x) + 2 x 2 ln(x)) 

3(l-x) 4 ' 
x(22 - 20x - 2x 2 + 16x ln(x) - x 2 ln(x) + 9 ln(x)) 



3(1 -x) 



with the corresponding mass ratios 
x t = mf/m 



x, 



2 / 2 



X 



99 



2 / 2 

m 5 /m.- . 



(19) 

(20) 
(21) 

(22) 



Bo(x t ) and Co(x t ) are the box and Z° penguin diagram functions in the Standard 
Model respectively. The function H Q (x qx ) appears in the SUSY contribution to the 
sdZ vertex [[RJ. The functions F (x gq ) and Go(x gq ) enter the contributions of 7- 
magnetic and chromomagnetic penguin operators respectively 



3.3 Effective couplings 

The SUSY Wilson coefficients which we have given above depend explicitly on the 
sparticle masses via the functions H , F and G . The dependence is not very strong, 
as can be seen from Fig. 0, where we plot the three functions normalized to their 
values at x = 1 (H (l) = -1/96, F (l) = 2/9, G (l) = -5/18). On the other 
hand the relations between e'/e and the rare decays which we want to investigate 
here, are almost independent from the spectra of the SUSY particles. In fact these 
relations are most conveniently described in terms of three effective couplings defined 
as follows: 



At 



(23) 



It is worthwhile to point out that most of the results presented in Section |6| are 
valid also if these couplings are defined in a more general way, starting from the 
Wilson coefficients of Z-penguin and chromomagnetic operators. This way one could 
efficiently include also subleading contributions in the mass-insertion approximation. 
This is however beyond the scope of the present analysis. 
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0.2 



1.8 



2.6 

x 



3.4 



4.2 



Figure 1: Dependence on x of the functions H (x)/H(l) (solid), G(x)/G(l) (dashed), 
F(x)/F(l) (dot-dashed). 

4 Basic Formulae for e' je and Rare Decays 

In this section we collect the formulae for e' je and rare K decays which we have 
used in our analysis. These formulae can be considered as the generalization of the 



corresponding expressions in [15] to include contributions of the chromomagnetic and 
7- magnetic operators to e' je and Kl — > 7r°e + e~ respectively. However, we stress that 



here we will treat the effective sdZ vertex differently than in [|15j , separating explicitly 
SM and supersymmetric contributions as shown in Q). The latter will be described 
in terms of the effective coupling A t defined in (^3|) . 



4.1 Magnetic contributions to e'/e and Kl — > ir°e + e 

The matrix elements of the magnetic operators between a K° and an n-pion 
state are difficult to calculate. In the following we will normalize them by using 
the value obtained in model calculations, and introduce the corresponding B factors 
which we will then vary inside our estimates of the uncertainties. We will use: 

<0nr)«l<SI*°> = f—^-ml B c , (24) 
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(AQ+\k°) 



Q d e iV2 
I07r rriK 



((nn) I=0 \Q+\K°) = (n°\Ql\K ) = 



(25) 
(26) 



For Bg = 1 Eq. (^4j) corresponds to the result of Ref. |2(J obtained at leading 
nontrivial order in the chiral quark model. We remark that the m 2 suppression 
of the matrix element is valid only at this order, and that terms proportional to 
m 2 K arise at the next order both in the 1/N C and in the chiral expansion. Large 
corrections to Bq = 1 are therefore rather plausible, and to take them into account 
we will use in what follows \B G \ = 1—4. As for Bt-, a value very close to one can be 
obtained for instance in the framework of vector meson dominance, as in |ZT| . Other 



estimates give very similar values (see e.g. ]22|]). As a conservative range of variation 
for this parameter we adopt \B?\ =0.5 — 2. Concerning the sign of Bt and Bq, the 
above model-dependent considerations indicate that it is positive in both cases. We 
stress, however, that this conclusion is not based on first principles. 
Using (|24"D we write the chromomagnetic contribution to e'/e as[] 

11V3 oj m^m 2 K a s {rrig) 

C\\ir - Yli'l 4,,i /•' liri rn , \ in 



Re It 



m~ g 



-r]B G ImA" 



(27) 



QAtt \e\ Re(A ) F w (m s + m d ) 

where rj contains the effect of the scaling from rrig down to m c (which is the scale 
at which the quark masses have to be given) and can be found in (|16|). Using 
a s (Mz) = 0.119 we then obtain 



Re f- 



G 



where 



R,, 



a s {m~ q 



209i? g ImA- 



23 

51 500GeV 



(2i 



R S B G . 



_a s (500GeV) 

As for the magnetic contribution to the direct CP-violating component of 
7r°e + e - , we notice that by using Eq. (p5|) one can write 



(29) 



(7r°e + e-|Q+|^°) 



Qd<yB T 



(7r°e + e-\Q 7V \K°) 



(30) 



where Qrv{A) — {sd)<v-A){^)v{A)- Employing the notations of |l9l and dropping for 
a moment the supersymmetric contribution to Z ds we get 

, 2" 



BR(K L 



7r°e + e- 



)dir = 6.3 ■ 10' 



(lm\ t y 7A ) 2 + (lm\ t y 7V + ImA+y 7 
"In our conventions Re A = 3.326 • 10" 4 and F n = 131 MeV. 



(31) 
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where f-y^viA) is the Wilson coefficients of Qrv(A) (the numerical values can be found 
in [|19j) and y 7 is defined by 

QdBr 



ImA+?/ 7 



Im 



F o( x gg) 

lGo(Xg q ) 



where R as = a s (m,g) / 'a s (500GeV) . 



+ 8 1- 1.13R as 21 



(32) 



4.2 Supersymmetric e'/e 

We decompose the SUSY contributions to e'/e as follows: 



(33) 



Z \ / G 

where the first term is the contribution from the supersymmetric effective sdZ vertex 
and the second is the contribution of the chromomagnetic penguin operator already 
discussed and given in (|28|). 
From ITHl we have 



Re 



where 



(9.- 


\.2-R s \rf\Bf' 2) ' 


R s = 


158MeV 


m s (m c ) + m d (m c )_ 



ImA 



t ■ 



(34) 



(35) 



(3/2) 

and Bg is the usual non-perturbative parameter describing the hadronic matrix 
element of the dominant electroweak penguin operator. Finally \r z \ is a calculable 
renormalization scheme independent parameter in the analytic formula for e'/e in 
p| which increases with af^(M z ) and in the range 0.116 < o^(M z ) < 0.122 
takes the values 



For R s we will use the range 



7.1 < \rf\ < 8.4. 



1 < R, < 2 



(36) 
(37) 



which is compatible with the most recent lattice and QCD sum rules calculations 
as reviewed in |J. Note that R s is defined as in which differs from |^] where 
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(8) 

158MeV has been replaced by 137MeV. Correspondingly the updated values of \r z \ 
given in || have been rescaled appropriately. We consider the ranges in ( flop and 
(|37|) as conservative. Finally we will use as in || 



0.6 < Bf /2) < 1.0. 



(3* 



Our treatment of all the other parameters which enter in the SM estimate of e'/e 
will be explained in Section 6. 



4.3 Rare Decays 



Following [15] we have 



BR(K + -> tt h 



VV) 



BR+ M + 1.55 • 10" 



2X Re (A t A^) + 2A C Re A t + \ A t 



, (39) 



where BRg M is the Standard Model contribution given by 



BRg M = 1.55 • 10" 



(X ImA 4 ) 2 + (X ReAi + A c 



(40) 



where 



A c = -(2.11 ±0.30) ■ 10" 4 (41) 

represents the internal charm contribution [[24[] and Xq = Co — ABq = 1.52 is the 
combination of penguin and box diagram functions in ([[7]) evaluated at fn t (m t ) = 166 
GeV. For an updated discussion about the SM estimate of the branching ratio we 



refer to [ p5|j . 

Next, following [15] and including the contribution of the 7- magnetic penguin to 
Kt, — > 7i°e + e~ we have 



BR(K L -> ir°e + e-) 



dir 



BRg M + 6.78 ■ 10" 4 [2Xo Im X t Im A t + (Im A t ) 2 ] , (42) 



BR^ e M + 1.19 ■ 10" 



BR{K L 



A* A* Jsd 



2Y lm\ t ImA 4 + (ImA t ) 2 
+2.13 Im Ai(0.08 Im A t + 0.23 Im A+y 7 ) 
+ (0.08ImA t + 0.23ImA+y 7 ) 2 , (43) 

BR^ + 6.32- 10^ 3 [2(F ReA t + A c )ReA t 

+ (ReA t ) 2 }, (44) 
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where the Standard Model contributions are given as follows 

2 



BRg M 

-DJ^SM 



6.78- 10 



-4 



X Im X t 



1.19 • 10- 4 (ImA t ) 2 [^ 2 + (1-0 + 0.08C ) 
6.32- lQ~ 3 \Y Re\ t + \ l2 



Here Y = C - B = 0.97, C = 0.79 and 



A c = -(6.54 ±0.60) ■ 10" 



(45) 
(46) 
(47) 



(48) 



represents the charm contribution to Kl — > [24 . 

Using (|39"D, ( f42| ) and (|44] ) one finds the following useful formula |15 



BR^ 



7T VV\ 



1.55-10" 



±3.97x/k • 10" 4 - 3B Re X t + A c 



+0.229 • BR(K L ttW 



(49) 



where 

A c = A c - A c 
and k is defined through 



;i.46±0.30) • 10" 



BR(K L -> y u + /i-) SD = « • 10" 



(50) 



(51) 



In evaluating A c we have included the correlation between A c and A c due to their 
simultaneous dependence on and m c P^| . The upper bound on BR(K + — > 
-n + vi>) is obtained for negative sign in (|49|) which corresponds to ReA 4 < C \ Re A t | 
(or ReZ sd < 0). 



5 Indirect bounds on supersymmetric contribu- 
tions 

5.1 Preliminaries 

We now discuss the presently available constraints, not directly obtained by e'/e 
or rare decays, on the flavour- changing mass insertions introduced in Section [3]. 
A general model-independent constraint on left-right mass insertions is dictated by 
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vacuum stability. In particular, the requirement of avoiding charge- or color-breaking 



minima or unbounded- from-below directions in the SUSY potential implies |26| 



12(21) 



\/3' 



m.. 



rrig 



(3 



V3m t 



(52) 



Given the large difference between top and strange quark masses, the two constraints 
in ( [52]) are numerically very different. However, when translated in bounds for the 
corresponding contributions to e'/e they look rather similar. Neglecting the depen- 
dence on the sparticles mass ratios, that is rather mild, we obtain 







|A t | 


< 10- 500GeV . 









A* I < 3 • 1CT 3 



/500GeVV 



(53) 



which leave open the possibility of large contributions to e'/e (up to ~ 10 2 ) both 
from ImA~ and ImA t . Concerning the bound on ImA+ relevant to Kl 



7r°e + e 



we further note that up to unlikely cancellations among (5^) 12 and (S^ R ) 2 i one 
expects 



ImA" 



ImA; 



Indirect bounds on and At can also be obtained by |AS| 



(54) 

2 amplitudes, 



barring the possibility of accidental cancellations. In the case of {5^ R 



12(21), 



the 



indirect constraints imposed by ex and Am^ are rather mild uM and essentially 



do not affect the bound in (|52j). In the case of A t , the constraints from \AS\ = 2 
amplitudes are of two types: those imposed by chargino box diagrams ]16[[| and those 
obtained via radiative corrections, relating (5^)23(^^)13 t° (^ll)i2- It turns out that 
the constraints via radiative corrections using Renormalization Group evolution are 
more severe than the ones from chargino box diagrams. We therefore discuss the 
former constraints in some detail. 



5.2 Bounds on A f via Renormalization Group 

The presence of a large double mass- insertion of the type (uj j — x — u s L ) could 

have a sizable indirect effect on the mixing of the first two generations, that is strongly 

t We note that the chargino contribution to \AS\ = 2 amplitudes has been overestimated in |lq| 
due to a missing factor 1/4 in the r.h.s. of Eq. (3.4). Moreover, though formally correct, Eq. (3.5) 
of [jl6| does not correspond to the expansion of Tt\AS\=2 near Xki = 1 (due to the missing factor 
1/M~ k ). Taking into account these two corrections, we found that the bounds on At in Eqs. (3.6-7) 
of pffl should be increased by a factor 3. 
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<v> <v> 

X X 



Figure 2: Diagrams through which the trilinear A u couplings may generate a 
Sl(r) -> <Il{r) transition. 



constrained in the down sector [13]. Indeed, the trilinear couplings A"' induce at 
one loop a flavour- changing mass term for both left- and right-handed squarks, i.e. 
give a radiative contribution to the off-diagonal elements of the mass matrices nig, 
m| and m| [2~7||. The diagram which generates such an effect is depicted in Fig. 



0, together with the diagram with the double LR mass-insertion which yields the 
d a — sa (A = L, R) transition. A naive order-of-magnitude comparison between the 
two diagrams (say, at low momentum q 2 flowing along the squark line) would lead 
one to say that the loop diagram is suppressed with respect to the tree one by a 
factor Mg/(16ir 2 (v ) 2 ) ~ 10~ 2 , if we assume that Ms is not much bigger than the 
electroweak-breaking scale. However, this suppression factor, which dominates over 
the finite part of the loop diagram, can be balanced by a large logarithm arising in 
the divergent part of the diagram. In particular, in a scenario with Mx ~ 10 16 GeV, 
the loop diagram yields a large logarithm of the form ln(M x /Mg) ~ 64 for Ms ~ 10 2 
GeV, therefore compensating almost completely the suppression factor. 

To bring this discussion on more solid grounds the tools to be used are the 
renormalization group equations (RGE) for the soft SUSY-breaking couplings p8| . 
If we neglect all entries in the Yukawa matrices but y t and y&, the RGE for the (1,2) 
matrix element of nig reads as follows (t = In M x /q 2 ): 



(55) 



where the ellipsis stand for terms which, according to the vacuum stability bounds 
are suppressed by (mb/m t ) 2 at least. Let us now imagine for a moment that the A" 
matrix elements do not evolve. Then we get for the radiatively generated part of 
(m Q )i 2 : 

( m Q)l2 i M S) = A 13 A 23 > (56) 
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that, when translated into the usual 5's becomes (for Ms = 300 GeV and Mx 
2 • 10 16 GeV, and tan/3 = 5): 



(CDi2 d = 1-3 • {5 u LR hM R )lz ■ (57) 

(A similar expression can be obtained for the couplings). If both the LR 

couplings were close to the vacuum stability bounds, this contribution would be of 
order one, and would violate the bounds which were obtained by comparison to the 



phenomenology of the AS = 2 transitions JTj]. By reversing the argument, and 
assuming there is no cancellation with the initial value of (nig) 12 at Mx we can 
obtain a bound on the product of the two LR couplings. 

In order to obtain the correct numerical value for this bound we have to do 
a complete calculation and take into account also the evolution of the A" matrix 
elements. In the same approximation as above (i.e. keeping only the y t and ?/& 
entries in the Yukawa matrices, and neglecting all the A u,d matrix elements whose 
vacuum-stability bound is not proportional to m t ), the RGE for the A" matrix 
elements read as follows: 



dt 8tt 



16 13 T 

—a 3 (t) + 3a 2 (t) + — ax(t) - -^W)? 



A£ (i^3) . (5* 



The one-loop evolution of the Yukawa coupling and of the gauge coupling constants 
in the MSSM is well-known, and can be found, e.g., in [p8[| . The boundary conditions 
which we have used for these equations are the following (for the scales Ms = 300 
GeV and M x = 2 ■ 10 16 GeV): 

y t {M s ) = 0.92 ±0.03 , 
y b (M s ) = 0.084 , 

ai (Mx) = 0.040 ±0.001 (i = l,...,3). (59) 

For simplicity we have evolved back from Mx all three gauge couplings from their 
unification value. With these boundary conditions, the solution of the RGE equation 
for (mg) 12 (Ms) is the following: 

(m|) 12 (M s ) = (m Q ) 12 (M x ) - K HM ^^ ] (A^A^) (M s ) (60) 
K = (0.67 ±0.05) , 
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where the uncertainty is mainly due to the top mass. As it is seen, the simplified 
solution (|56|) is numerically not very different from the complete one in (|60|). It 
is interesting to note that also here the large top mass plays an important role: 
the Yukawa coupling largely compensates the effect of the gauge couplings in the 
evolution of the A" 3 matrix elements. Neglecting the Yukawa term in (|58|) , the 
numerical coefficient —0.67 goes down to —0.34. Disregarding the unlikely possibility 
of a strong cancellation between the two terms on the r.h.s. of flBTf ) we can obtain 
a bound for A t (for the numerical estimate we use again tan/3 = 5 and Ms = 300 
GeV): 

16n 2 sm 2 (3 



Ilm A f 



< 



\Hn(X n 



K\n{M 2 x /M 2 s )M s 



mm 



in ||lm((5 



ll) 12 



Im(5 



u \ 

LL)12 



(1.2 ±0.1' 



Ho(x qx 



lm(5 



ll) 12 



(61) 



and analogously for the real part. 

The left-left mixing among the first two generations of down-type squarks is 
strongly constrained since it appears in gluino- mediated |AS| = 2 amplitudes fT4"] . 
Since (S^ L ) \i enters quadratically in |A£| = 2 transitions, one gets the following 
bounds from AMk and Ek respectively fL4| : 



|Re(5 



LL)\2\ 



|Im(tf 



LL) 12 I 



< 2.4-10" 



< 1.9-10" 



\ 



4/ 6 (l) + ll/ 6 (l) 



4:Xg q f 6 (Xgq) + 11/ 6 (X 



<}<-}) 



\ 



4/ 6 (l) + ll/ 6 (l) 



^ gq f 6 (Xgg) + Uf 6 (x 



300GeV 



300GeV 



(62) 
(63) 



where the functions / 6 and / 6 are defined in [|14|. The combination 4x/ 6 (x) + llf e (x) 
has a zero at x = 2.43, so that close to this particular value of the gluino-squark mass 
ratio the bounds (|62] - |63D become irrelevant. On the other hand, this value is excluded 
in the present scenario where Mx ~ 10 16 GeV, because the evolution of the masses 
via RGE down to electroweak scales gives the condition x gq < 1.3 for the scalars of 
the first two families [^]. Moreover, even if the limits coming from gluino exchange 
could be evaded, the analogous limits coming from chargino exchange, which are not 
much weaker, would still hold. 

Using Eqs. (|6T|-|63D it is possible to obtain bounds on lmA t that are more strin- 
gent than the one in Eq. fl53|). However, the precise size of these constraints depends 
strongly on the phase of A t : if the double insertion is purely imaginary, the constraint 
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from Ek is ineffective and Im A f can be substantially larger than in the case in which 
ReA t is different from zero. 



5.3 Scanning of the SUSY parameter space and model-de- 
pendent considerations 

Taking into account the analytic bounds discussed so far, we will now proceed esti- 
mating the maximal allowed size of Im A t in terms of various SUSY parameters. To 
do so, one has to face the usual problem of scanning efficiently the parameter space. 
In this particular case, the phases of the relevant FCNC parameters are crucial: as 
we stressed above, the stringent constraint from ek is not effective on pure imaginary 
(double) mass insertions. 

To obtain an estimate of model-independent limits on SUSY contributions, we 
scan randomly with uniform distribution the parameter space corresponding to a 
reasonably natural determination of Mz- More precisely, we choose the relevant 
parameters in the following intervals: —300 GeV < fx < 300 GeV[j], 100 GeV < 
M 2 < 250 GeV, 3M 2 < m Qi2 < 5M 2 , M 2 < m Li2 < 2M 2 , 0.4 < m- tR < m Qi% . 
Moreover we assume unification of gaugino masses and we discard points in which 
(M 3 /rriQ i2 ) 2 > 1.3, the charginos are lighter than 90 GeV, the charged sleptons lighter 
than 80 GeV or the gluinos lighter than 180 GeV. The limits we get however do not 
significantly depend on the details of the scanning procedure. We focus here only 
on the possibility of large enhancements with respect to the SM due to the double 
mass insertion contribution to Im Z^. Since the effects of single mass insertions have 



already been analyzed in detail in Ref. |18| and have been shown to be smaller, or 
at most of the same size of the SM contribution, we do not take them into account 
in the present analysis. 

As we discussed before, the most stringent upper limits on the double mass 
insertion come from Ek and Am^ through the RGE evolution. To estimate the 
maximal possible effects, we first choose the double mass insertion phase, then we 
choose the corresponding absolute value as high as the highest limit found scanning 
the parameter space. In Figure |3|, we plot the maximal possible value of | ImA t | as 
a function of argA t . It is evident that the stringent constraint from Ek forces ImA t 



■^We use a real fi to avoid problems with the electric dipole moment of the neutron. 
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Figure 3: Limit on | ImA^ imposed by Am^ and Ek, through RGE evolution, as 
a function of argA t . The dashed line shows the SM contribution to ImZ^s for 
ImAt = 1.33 ■ 10" 4 . 

to be smaller than or of the same order of the SM contribution to lmZ^ s , unless 
argA t is very close to ±7r/2. Therefore a large enhancement of Im^ s with respect 
to the SM can only happen if the double mass insertion is large and almost purely 
imaginary. In this particular case, combining flST]) and fl62|) we can write 



ImAJ < 3 • 1(T 4 



H a (x. 



\ 



4/ 6 (l) + 11/ 6 (1) 



300GeV 



rrig 



4:X gq f 6 (x gq ) + Uf e (x gq ) 

As we shall discuss in the next section, this particular case can be tested experi- 
mentally in a clear way by studying rare K decays: if for example BR(i^£ — > -n^vv) 
will be found to agree with the SM expectations, then the possibility of a large Im A t 
will be ruled out. 

The constraints we considered on the relevant mass insertions can be evaded in 
corners of parameter space, but this holds only if an unlikely fine-tuning is allowed. 
For example the limits from Amx and sk can be evaded if there is a cancellation 
among the different supersymmetric contributions to them, or the limit from RGE 
can be evaded if there is a cancellation between the initial value of the insertion and 
the RGE contribution. 

Since the insertions are pushed up to their experimental limits the results plotted 
should not of course be considered as predictions but just as maximal possible effects. 
Our framework is in fact general enough to include any supersymmetric extensions 
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of the SM with minimal field content. This on one hand insures that we are not 
missing potentially large effects. On the other hand, one might ask whether values 
of |ImA t | as large as those ones in the shaded region of Fig. |^ naturally arise in 
super symmetric models. Unfortunately, within the most common models this is not 
the case, as we will now briefly show. 

Explicit models account for the strong constraints on soft supersymmetry break- 
ing terms in different ways. In some cases the mechanism communicating the su- 
persymmetry breaking guarantees that FCNC and CP-violating processes are under 
control. This is the case e.g. of gauge mediated supersymmetry breaking and of 
minimal supergravity (SUGRA). In other cases, further ingredients are necessary. 

In the minimal situations, a quick estimate yields 

where Mx ~ 10 16 has been used to estimate A^ and m' s , m'g are dimensionful 
combinations of diagonal soft parameters. Eqs. (|65D and (^6|) show that A t and A^ 
give rise to negligible effects compared to the SM ones. 

On the other hand the universality hypothesis used in minimal SUGRA has 
not a compelling justification. In this and other cases in which the mechanism 
generating the soft terms does not guarantee that FCNC are under control, the 
potential FCNC problem must be solved by further symmetries. From this point 
of view the issue of why the scalar mass eigenstates are so degenerate or so aligned 
with the corresponding fermion eigenstates is the supersymmetric version of the 
issue of explaining the structure of fermion masses and mixings. If the latter is 
accounted for by flavour symmetries acting on the fermion generation indices, in a 
supersymmetric theory the same symmetry acts on the corresponding scalar indices. 
As a consequence, whatever is the symmetry, since the Yukawa and the corresponding 
soft trilinear interactions have the same quantum numbers, the structure of their 
coupling matrices is the same. Within this class of models it is therefore possible 
to show that the LR mass insertions involving the third generation, and in turn the 
double insertion, are similar to those obtained in the minimal models. This is not 
what happens for (<5f/j)i2, that can be shown to be of the right order of magnitude 
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to generate the experimental value of e'/e [11]. Therefore the most likely situation, 
as far as the most common SUSY models are concerned, is somewhere between the 
case of A t ~ and A g ^ and the case of A t = A g = 0. We stress, however, that the 
flavor structure of the supersymmetry breaking is far from having been established. 
It is then worthwhile to investigate also more exotic possibilities, like the one of a 
large ImA 4 , as far as these are not ruled out by phenomenological constraints. 



6 Numerical Analysis 

6.1 Strategy 

We are now ready to discuss magnitude and relations among possible supersymmetric 
contributions to e'/e and rare decays. To this purpose it is useful to distinguish 
between three basic scenarios: 

Scenario A: [ImA t = 0, ImA^ ^ 0]. 

This scenario is close to what happens in most SUSY models since, as we have 
seen in the previous section, the sdZ vertex can receive sizable corrections only 
in a specific region of the parameter space. In this case e'/e can be affected 
only by the chromomagnetic operator and, as shown in Section [PL among the 



rare modes only Kl — > n e + e~ is sensitive to this SUSY contribution. On the 
other hand if ReA t is substantially different from zero also K + — > 7i + uu can 
be significantly affected. 

Scenario B: [lmA t ^ 0, ImA^ = 0]. 

In this scenario the possibility of large corrections to e'/ e is not favoured from 
the point of view of the parameter space, but is an interesting possibility to 
be investigated in a model-independent approach. If this is the case, sizable 
effects are then expected both in Kl — > vr°z/z/ and Kl — > 7r°e + e~. 

Scenario C: [ImA t ^ 0, ImAj ^ 0]. 

This represents the most general case. Note, however, that the requirement 
of having sizable cancellations in e'/e, between supersymmetric contributions 
generated by the chromomagnetic operator and the sdZ vertex, implies an 
additional fine-tuning with respect to scenarios A and B. 
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We will also follow |T5) and consider three scenarios for \ t , which enter Standard 
Model contributions and its interference with supersymmetric contributions to rare 
decays and e'/e. Indeed there is the possibility that the value of X t is modified by 
new contributions to e and B^ s — B^ s mixings. We consider therefore three scenarios: 

• Scenario I: At is taken from the standard analysis of the unitarity triangle 
and varied in the ranges: 

1.05 • 10~ 4 < Im X t < 1.61 • 1(T 4 (67) 

2.3- 1(T 4 < -Re At < 3.8- 1(T 4 (68) 

• Scenario II: Im X t = and Re At is varied in the full range consistent with 
the unitarity of the CKM matrix: 

1.61- HT 4 < -Re A t < 5.6- 10" 4 (69) 

In this scenario CP violation comes entirely from new physics contributions. 

• Scenario III: Xt is varied in the full range consistent with the unitarity of the 
CKM matrix: 

- 1.73 • 1(T 4 < ImAi < 1.73 -KT 4 (70) 
This means in particular that Im At can be negative. 

We would like to emphasize that the scenarios II and in particular III are very 
unlikely and are presented here only for completeness. We stress that if one uses 
the Standard Model expressions for B° — B° mixings, Ek and sin 2(3 one gets results 
for the CKM matrix which are compatible with the |Vub/V^| constraint, which is 
insensitive to new physics. In view of the coherence of the Standard Model picture, 
corrections to the processes in question so large as to make Im X t negative, or Re A t 
way outside the range in Eq. ( |B8"D look rather improbable. We believe that if the new 
physics has an impact on the usual determination of At, the most likely situation is 
between scenarios I and II. 



6.2 e'/e 

We shall now proceed extracting ranges for the effective SUSY couplings from the 
experimental data on e'/e in the basic scenarios A-C defined above. These will then 
be used to estimate the branching ratios of the rare decay modes. 
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Figure 4: Linear relation between A g and At for Re(e'/e) SUSY = 2 • 10 3 . The solid 
line is for {Bf /2 \ R s , rf} = {0.8,1.5,7.8}, the dot-dashed for {sf /2) , R„ rf} = 
{1.0, 2.0, 8.4} and the dashed for {B^ /2) , R s , rf} = {0.6, 1.0, 7.1}. The vertical lines 
show the RGE bound flBlf) for rriq = 300 GeV and {x qx ,x gq } = {3, 1} (dotted) or 
{x gx ,x gq } = {9,1.3} (dashed). 

Assuming that the SM contribution to Re(e'/e) is around its central value, as 
given in ||, and therefore much smaller than the experimental result, there is a lot 
of room for SUSY to contribute to this quantity. Detailed bounds on Re(e'/e) SVSY 
depend on the various parameters entering Re(e'/e) SM , as well as on the experimental 
result in (||), however, as a simplified starting point for our discussion, we assume at 
first 

(, \ SUSY 
j J =2-10~ 3 . (71) 

This value has to be taken only as a reference figure: it could be interpreted either 
as the difference between the experimental result and the SM contribution or as the 
true value of ~Re(e'/e) in the limit of a real CKM matrix. 

Since our formula for the SUSY contribution Eq. (^) contains only two free 
parameters, ImA 4 and ImA~, Eq. (|7T| ) defines a straight line in the (ImAt,ImA g ) 
plane, which represents the general solution within scenario C. This is shown in 
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Fig. |] for three different sets of {B^ 2 \ R SJ r^ 8 }. Decreasing the reference value in 
(ffl|) corresponds to a translation of the straight lines toward the origin; the intercepts 
of the lines with vertical and horizontal axes define the solutions within scenarios A 
and B, respectively. 

As it can be noticed, if A~ = 0, then lmA t must be negative, i.e. the SUSY 
contribution to the sdZ vertex must be opposite to the SM one in order to produce 
a positive contribution to e'/e. The minimum value of |ImA f | with A~ = is 
found for the maximum values of B^ 2 \ R s and r% . In this case SUSY and SM 
contributions to the sdZ vertex cancel almost completely and the experimental value 
for e'/e is roughly reproduced by QCD penguin contributions. On the other hand, 
the maximum allowed value of I ImA+l with A~ = is found for the minimum set 

I H g 

of {B^ 2 \ R s , Tz }. In this case the sdZ vertex has an opposite sign with respect to 
the SM case and is largely dominated by SUSY contributions > 6). This 

solution is still allowed by the RGE constraint (0), provided the sparticle masses 
are not too high. The situation of course changes if one allows also ImA ff to be 
different from zero. In particular, for large 10~ 5 ) and positive values of i? 9 ImA~ 
a positive ImA 4 is needed in order to avoid too large effects in e'/e. 

In the limit where the standard determination of the CKM matrix is valid, a 
quantitative estimate of the ranges for ImA~ and ImA 4 , within scenarios A and B, 
can be obtained by subtracting the SM contribution from the experimental value in 
(0). Following [Tj|, we parametrize the SM result for e'/e using the approximate 
formula 



e' 



Re I -J = Im A t 
with U 



■1A + R S 



l.l|4 8) |^ 1/2) + (1.0 - 0.67|4 8) |)^ 3/2) 



(72) 



Im X t = (1.33 ±0.14) ■ 1(T 4 . (73) 

Varying lm\ t and the experimental value @ within 2a intervals, choosing B^ 2 \ 

R s and r% as discussed in Section [4.2| and, finally, assuming 0.7 < Bi 1/2) < 1.3, we 
find: 

- f i \ SUSY , 

- 15.5 • < Re (M <30.1-10" 4 , (74) 

-9.3-10~ 4 < ImA t < 1.7 ■ 10~ 4 (ImA;=0), (75) 

-0.7-10~ 5 < .RglmA^ <1.4-10" 5 (ImA t = 0) . (76) 
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It is interesting to note that the range of ImA^ is well within the bound (|52|). 
therefore e'/e provides the most stringent bound on | ImA~| within scenario A. Sim- 
ilarly, e'/e provides the most stringent model-independent upper bound on ImAj 
within scenario B. On the other hand, the lower bound on lmA t imposed by e'/e is 
weaker than the bound (|64] ) for m q > 200 GeV and x gq < 1.3. 

To show the possible improvement due to more precise measurement of e'/e we 
show how ©-(0) are modified if we fix Re(e ; /£) exp = 20 • 10" 4 . We find 

. / / \ SUSY . 

-7.5-10" 4 < Re(|) < 19.7 • 10" 4 , (77) 
-5.9-10" 4 < ImA t <0.8-10~ 4 (ImA: = 0) , (78) 
-0.4-10~ 5 < R g lmA~ <0.9-10" 5 (ImA t = 0). (79) 

6.3 Rare Decays 

The rare decays Kl — > 7r°z/z/ and Kl — > 7r°e + e~ provide in principle a powerful tool to 
clearly establish possible SUSY contributions in CP- violating \AS\ = 1 amplitudes, 
and also to distinguish among the three scenarios introduced in Section [67 

6.3.1 Scenario A 

Within scenario A only Kl — > vr°e + e _ among these two modes is affected by SUSY 
corrections. Setting R as = 1 in (^) we can write 



Im A+y 7 = 35.5 R g Im A^ 



Im A+ 
Im Ajj 



0) 



where the numerical coefficient has been obtained for x gq = 1 and can increase at 
most to 37.0 if we impose x gq < 1.3. Assuming R g lmh~ = 10~ 5 , as obtained from 
Fig. |], and fixing to unit the two ratios among square brackets in (|80|), we obtain 
ImAj ji 7 = 3.5- 10~ 4 . Using this figure in (pE3"| ) we find that the additional contribution 
to BR(K L — > vr°e + e~)dir is positive and ranges between 3 and 4 in units of 10~ 12 , 
depending on the value of ImA^. This effect, which represents the typical size of 
the SUSY contribution to Kl 7r°e + e~ expected within scenario A, is certainly 
difficult to be observed. However, we stress that this conclusion depends strongly on 
the assumptions made for ImA+/ImA~ and B T /(B G ^/Rs). 

According to the ranges of Bt, Bq and R s discussed in Section |], we expect 

0.09 < — ^= < 2 . (81) 
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On the other hand, it is more difficult to estimate ImA+/ImA~ without specific 
assumptions on the SUSY soft-breaking terms. In minimal models it is natural to 
assume (5£r)i2 ^> (^lr)21j that implies 

ImA+ / X 

ImA ff v ' 

but we cannot exclude sizable deviations from this figure in generic scenarios. 

In Table [l] we report the upper bounds on BR(Kl — > 7r°e + e~)di r , for different 
values of the two ratios. To this end we have used the expressions for e'/e and 
BK(Kl — > 7r°e + e _ )dir given in Section 4 with the Standard Model contribution for 
e'/e given in ([72]). Scanning the parameters B^ 2 \ B^ 2 \ R s and rjj as discussed 
in Section [O] and 6.2, varying lmX t according to ( p7| ) (scenario I), we find the 
results in the third and fourth column which correspond to two choices of e'/e. As it 
can be noticed, results in the ball park of 1 cannot be excluded even under the 
assumption ([82]). 

The dependence of BR(A^ — > 7r°e + e~)dir on the value of Im At is shown in Table ^|. 
If the CKM matrix is real and |(Im A+/ Im Aj)B T /(B Gy /R^)\ > 1, we find BR(K L -» 
7r°e + e~)dir ~ few x 10 -12 , similarly to the SM case. On the other hand values 
substantially larger than 10 -11 are obtained within scenario III. Note, however, that 
the large results quoted for ImA^/ImA" < are very unlikely, since are obtained 
for the maximum negative value of ImAj. 

Concerning Kl — > n^vv, its branching ratio in scenario A stays close to the 
Standard Model value provided the usual determination of Im X t is not substantially 
decreased through supersymmetric contributions to Ek- Because of the unitarity of 
the CKM matrix Im X t can only be marginally increased over its SM value. On the 
other hand if Im A t = a clear signature for scenario A would be a vanishingly small 
BR{K L -> 7f°uu) (< I(T 14 |3). 

The case of K + — > n + uu is different as it is dominantly governed by Re At and 
Re Aj. The upper bound on BR(A' + — > 7i + uu) can be obtained by using equation 
(H) together with the bound g§ ^ f§ 

BR(K L -> m + /Osd < 2 - 8 ■ 10_9 (83) 

i.e. k = 2.8. Choosing then (— ReAt) max = 3.8 • 10~ 4 (scenario I for At), as obtained 
in the Standard Model, or (— ReA f ) max = 5.6 ■ 10 -4 (scenarios II and III), we find 
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Im A+ / Im A g 


Bt/ {B G \TR^) 


BR(K L -> vr°e + e-) dir 


BR(K L -> 7r°e + e-) dir 


-1 


1.0 


9.4- 10~ 12 


7.8- 10~ 12 


-1 


0.5 


7.8- 10- 12 


7.0- 10- 12 


-1 


1.5 


1.1 • 10~ n 


8.5- 10~ 12 


-2 


1.5 


1.8 • 10~ n 


1.1 • lO" 11 


1 


1.0 


1.3- lO^ 11 


1.0- 10" 11 


1 


0.5 


9.3- 10~ 12 


8.2 • 10~ 12 


1 


1.5 


1.8- 10 -11 


1.3- 10 -11 


2 


1.5 


3.7- 10~ n 


2.3 ■ 10~ n 



Table 1: Upper bounds on BR(K~l — > 7r°e + e )dir within scenario A, for different 
values of ImA+/ImAj and B T /(B G ^fR^) consistent with 12 < 10 4 Re(e7e) < 30.4 
(third column) and Re(e'/e) = 20.0 • 10 -4 (fourth column). The bounds are obtained 
setting x gq = 1.3 in order to maximize the numerical coefficient in (j30|). To maximize 
the interference of SM and SUSY amplitudes, R g Im A~ is chosen as the maximum 
(minimum) value allowed by e'/e for positive (negative) ImA+/ImA~. 



Im A+/ImA 9 


BR(K L -> vr°e + e-) dir (II) 


BR{K L -> vr e+e-) dir (III) 


-1 


1.8 (0.8) -10- 12 


2.5 (2.1) -10" 11 


-2 


7.3 (3.2) -10" 12 


5.7 (4.5) -lO" 11 


1 


1.8 (0.8) TO" 12 


1.5 (1.2) -10" 11 


2 


7.3 (3.2) -10- 12 


2.5 (1.7) -10" 11 



Table 2: Upper bounds on BR(K L — > 7r°e + e within scenario A for lm\ t = (II) 
and | ImAfl < 1.73 ■ 10 -4 (III). The bounds are obtained setting Bt/ '(BgVRs) — 1 
and imposing Re(e'/e) < 30.4(20.0) • 10" 4 . 
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respectively 



BR(A+ -> tt+z/z/) < 1.70 • lO" 10 + 0.229BR(fsT L -> vr°z/z/) , (84) 

BR(K + -> vr+z/z/) < 2.03 • 10~ 10 + 0.229BR(K L -> 7r°z/p) (85) 

As the second terms on the r.h.s of these bounds are very small in this scenario we 
find BR(K + -> vr+z/z/) < 1.7- 10 -10 and BR(K+ -> vr+z/z/) < 2.1 -HT 10 . These results 
are also obtained if Re A t is varied in the full range consistent with the bound ( j33|) 
and with the RGE constraint fl88|) with ImA t = 0. Evidently as (|8~4"D and ( jS5f) have 
been obtained without the constraint fl88|), what matters in this scenario is (|83|). 



6.3.2 Scenario B 

Being strongly sensitive to ImA 4 and insensitive to ImA^, Kl — > is^vv represents 
the golden mode to identify scenarios B and C. We first discuss scenario B which 



corresponds to the case analyzed in |15| . This time, however, the effective sdZ vertex 
is additionally constrained by the renormalization group analysis of Section 5. 

The dependence of BR(A^ — > -n^vv) on ImA t is shown in the left plot of Fig. [|. 
As can be noticed, large enhancements with respect to the SM case are possible, but 
on the other hand we cannot exclude a destructive interference among SUSY and 
SM contributions leading to strong suppression of BK(Kl — > tx°ui>). 

If the standard determination of Im X t is valid, Eq. ( fj2p implies that BR(i^ — > 
7r°z/i/) can be enhanced with respect to the SM case only if 

ImA t < -2X ImAt or ImA t > . (86) 

The second possibility is excluded within scenario B if we require a positive SUSY 
contribution to e'/e. This is clearly shown by the second plot in Fig. which 
illustrates the relation between BR(A^ — > 7r°i/z/) and the SUSY contribution to 
e'/e within scenario B, assuming the standard determination of lm\ t . In this case 
large enhancements of BK(Kl — ► tt°uu) are possible, but only if R s and B 8 are 
close to their minimum values. On the other hand, if R s and/or B$ are large, then 
BK(Kl — > tx q vv) is more likely to be suppressed rather than enhanced with respect 
to the SM case. 

In order to be more quantitative we consider the three scenarios for X t defined 
at the beginning of this section. Next, as discussed in Section 6.2, ImAj can be 
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Figure 5: BK(Kl — > it^vv) as a function of Im A t (left) or as a function of (e' ' /s) z 
(right). In the left plot the solid (dot-dashed) parabola is for ImA 4 = 1.33 ■ 10~ 4 (0) 
and the vertical lines show the RGE bounds as in Fig. |j. In the right plot the three 
parabola are for ImA t = 1.33 • 1(T 4 and {B§ /2) , R s , rf} = {0.6,1.0,7.1} (solid), 
{B$ m ,R a ,rf} = {0.7,1.0,7.8} (dot-dashed) or {Bf /2) ,R s ,rf} = {0.8,1.5,7.8} 
(dashed). In both cases the horizontal lines denote the SM range of BR(Kl — > 7r°z/P) 
for 1.05 < 10 4 ImA 4 < 1.61. 

best bounded by e'/e and the renormalization group analysis of Section 5. ReA t 
can be bounded by the present information on the short distance contribution to 
Kl —>■ fi + fi~ and also by the RG analysis of Section 5, as we will state more explicitly 
below. These bounds imply a bound on BK(K + — > -k + vv). Since BR(fT + — > 7i + uu) 
depends on both Re At and Im At also the bound on Im At matters in cases where it 
is substantially larger than the Standard Model contribution to Im Z^. 

The branching ratios BK(Kl — > 7c°i>u) and BR(K^ — > 7r e + e _ ) ( j ir are dominated 
by (ImZ S( j) 2 . Yet, the outcome of this analysis depends sensitively on the sign of 
Im Z s d- Indeed, Im Z s d > results in the suppression of e'/ e and since in the Standard 
Model the value for e'/e is generally below the data, substantial enhancements of 
ImZ S( j with ImZ S( j > are not possible. The situation changes if new physics 
reverses the sign of ImZ S( j so that it becomes negative. Then the upper bound 
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on — ImZgd is governed by the upper bound on e'/e and with suitable choice of 
hadronic parameters and Im X t (in particular in scenario III) large enhancements of 
— Im Z S( i and of rare decay branching ratios are in principle possible. The largest 
branching ratios are found when the neutral meson mixing is dominated by new 
physics contributions which force ImAi to be as negative as possible within the 
unitarity of the CKM matrix. As we argued above, this possibility is quite remote. 
However, if this situation could be realized in some exotic model, then the branching 



ratios in question could be very high as demonstrated in [15]. 

In this context it is interesting to observe that in the case of supersymmetry 
such large enhancements of — Im Z S d while allowed by e'/e are ruled out by the 
renormalization group bound on ImAf considered in Section 5. As we will see in a 
moment the imposition of the bound (see Fig. [3]) 

|ImA t | < 5.0- 10" 4 (87) 



has in the case of a negative ImA ( a very large impact on the analysis in sup- 
pressing considerably the upper bounds on rare decays obtained there. 

In Table |] we show the upper bounds on rare decays for lmA t > for three 
scenarios of ImAf in question and two different lower bounds on e'/e. To this end 
all parameters relevant for e'/e have been scanned in the ranges used in scenario 
A except that Im A^ have been set to zero. In Table |] the case Im A t < for two 
different upper bounds on e'/e is considered. In the last column we always give the 
upper bounds obtained in the Standard Model. 

The inspection of Table ^| shows that only moderate enhancements of branching 
ratios are allowed by e'/e if ImA ( > 0. Moreover the case lm\ t = is excluded by 
the positive value of e'/e. If Im A t < 0, substantial enhancements of BR(Kl — > -n^vv) 
and BR(i^i — > 7r°e + e~)di r are possible as seen in Table |j. In particular in scenario 
III both branching ratios can be enhanced by one order of magnitude over Standard 
Model expectations. On the other hand the imposition of the the RGE bound (|S7|) 
plays an important role in this analysis. In Table [| we show what one would find 
instead of Table [|, for Re(e' / e) mauX = 30.0 • 10~ 4 , if the bound ( ]8"T| ) had not been 
imposed. Table [5] corresponds to the analysis in jl5] and shows very clearly that 
without the bound (1871) very large enhancements of branching ratios in question are 



( 3 /2) 

possible. One should note the strong sensitivity of the results to the choice of Bg 
in Tables |3| and where the bounds are governed by e'/e. On the other hand this 
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10 4 Re(£V^)min 


Scenario 


for A t : 


I 


II 


Ill 


SM 




10 10 BR(K L - 




1.2 (0.6) 


— 


1.4 (0.8) 


0.4 


12.0 


10 11 BR{K L - 


-> 7T°e + e _ )dir 


1.7 (0.9) 


— 


2.1 (1.1) 


0.7 




10 10 BR(K + - 


-> 7T + Z/Z/)* 


2.0 (1.8) 


— 


2.4 (2.2) 


1.1 




10 10 BR(K+ - 


-> 7T + Z/Z/) 


1.7 (1.7) 


— 


2.1 (1.9) 


1.1 




10 10 BR(K L - 




0.7 (0.4) 




0.9 (0.5) 


0.4 


20.0 


10 11 BR(K L - 


-> 7T e + e _ )dir 


1.1 (0.7) 




1.3 (0.8) 


0.7 




10 10 BR(K+ - 


-> 7T+Z/Z/)* 


1.9 (1.8) 




2.2 (2.2) 


1.1 




10 10 BR(K + - 




1.7 (1.7) 




2.0 (1.9) 


1.1 



Table 3: Upper bounds for the branching ratios of the rare decays Kl — ► it vis, 
Kl — > 7r°e + e~ and K + — > 7i + uu in the case lmA t > 0, ImA^ = 0. The results 
have been obtained in various scenarios for A t by imposing Re(e'/e) > 12.0 • 10~ 4 or 
Re(e'/e) > 20.0-10" 4 , with B^ 3/2) = 0.6(1.0). The * means that the BR(K + -> vr+z/z/) 
has been calculated using the bounds fl8~4j) and (p5|). Otherwise, the more stringent 
bound due to RGE, Eq. (|88|) , has been used. 



10 4 Re(£7e) max 




Scenario 


for A t : 




I 




II 


III 


SM 




10 10 


BR(K L - 




0.8 


(0.8) 


1.7 


(1.7) 


4.0 (4.0) 


0.4 


30.4 


10 11 


BR(K L - 


-> 7T°e + e")dir 


2.0 


(2.0) 


3.0 


(3.0) 


5.9 (5.9) 


0.7 




10 10 


RR{K+ - 


^ TI + VV)* 


1.9 


(1.9) 


2.4 


(2.4) 


2.9 (2.9) 


1.1 




10 10 


BR(K + - 


-> 71 + UU) 


1.7 


(1.7) 


2.1 


(2.1) 


2.7 (2.7) 


1.1 




10 10 


BR(K L - 


-> 7T°b>v) 


0.8 


(0.4) 


1.7 


(0.8) 


4.0 (3.8) 


0.4 


20.0 


10 11 


BR(K L - 


■> 7r°e + e~)dir 


2.0 


(0.7) 


3.0 


(1.4) 


5.9 (5.7) 


0.7 




10 10 


BR(K + - 


-> 7T + Z/Z/)* 


1.9 


(1.8) 


2.4 


(2.2) 


2.9 (2.9) 


1.1 




10 10 


BR(K + - 


-> 7T + Z/Z/) 


1.7 


(1.7) 


2.1 


(1.9) 


2.7 (2.6) 


1.1 



Table 4: Upper bounds for the branching ratios of the rare decays Ki — > 7r°z/z/, 
Ki — > 7r°e + e _ and K + — > 7i + uu in the case ImA t < 0, ImA^ = 0. The results 
have been obtained in various scenarios for At by imposing Re(e'/e) < 30.4 • 10~ 4 or 
Re(eVe) < 20.0- 10" 4 , with Bf /2) = 0.6(1.0). For an explanation of the * see caption 
of Table |. 
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Scenario for X t : I II III 


SM 


10 10 BR(K L -> n°vis) 


3.9 


6.5 


17.6 


0.4 


10 11 BR(K L -> 7r°e + e-) dir 


7.9 


11.5 


28.0 


0.7 


10 10 BR(AT+ -> tt+z/P) 


2.6 


3.5 


6.1 


1.1 



Table 5: Upper bounds for the branching ratios of the rare decays Kl — * 7r°z/z/, 
Ki —>■ TT°e + e~ and K + — > n + uu in scenario B, without imposing the RGE constraint 
(|7|) and using /2) = 0.6. 

sensitivity is absent in Table £| for Re(e' /e) max = 30.0 ■ 10~ 4 and in scenario III for 
Re(e'/e) max = 20.0 ■ 10~ 4 , where the bounds on K L — > ix^vv and K L — > 7r°e + e _ are 
governed by the renormalization group bound fl8"7p. 

Next we should make a few remarks on K + — > 7r + z/z/. The bounds on BR(fC + — > 
7r + z/i/) denoted by "*" in Tables |3] and § have been obtained by using the bounds ( jS4|) 
and ( p5]) for scenario I and scenarios (II, III) respectively. It should be emphasized 
that these bounds are rather conservative as they take only into account the RGE 
bound in ImA t (through Kl — > -n^vv) and the bound on ReA 4 from (|83|). On the 
other hand, if A 4 is almost purely imaginary, as required by the RGE constraints for 
a large ImA t , the upper bound on ReA t is generally stronger than the one from ( j33|) 
and one has milder enhancements of BR(f^ + — > tt + vv) than in the "*" case. That 
is, in order to find the true bound, the correlation between Im At and Re A t through 
RGE should be taken into account. In order to investigate this correlation we have 
repeated the analysis for K + — > 7i + uu imposing instead of ( [37]) the more general 
RGE constraints 



|A t | < 5.0 • 10~ 4 , | Re A t ImAt | < 0.8 • 10 



-9 



derived from ( |B"T| - |6"3{ ) . The results of this analysis are represented by BK(K + — > 
7i + uu) without "*" in tables. As expected the bounds are stronger than previously 
obtained. Moreover the sensitivity to e'/e diminished and the bounds are mainly 
governed by K L — > yU + /i~ and RGE. 
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Scenario for A t : 


T 
1 


TT 
11 


TTT 
in 


Olvl 




-> 7T f v') 


O.O 1 O.O J 


i. f i^i. g 




n a 


10 10 BR(fT+ - 


-> 7T + Z/Z/)* 


9 ft (9 

Zr.U ^Zr.U J 


9 4 (0 A\ 




1 1 
1 . 1 


10 10 BR(AT+ - 


-> 7T + Z/Z/) 


1.8 (1.8) 


2.1 (2.1) 


2.7 (2.7) 


1.1 


10 11 BR(K L - 


-> 7r°e+e-) dir [+] 


10.0 (9.3) 


5.7 (5.3) 


10.3 (9.7) 


0.7 


10 11 BR(AT L - 


-vr°e + e-) dir [-] 


5.7 (5.5) 


4.9 (4.5) 


6.8 (6.1) 


0.7 



Table 6: Upper bounds for the branching ratios of the rare decays Kl — > ~K Q vv y 
Ki — > TT°e + e~ and K + — > 7r + z/z/ in scenario C, for ImAt > and -R^ImA^ > 0, 
imposing Re(e'/e) < 30.4(20.0) ■ 10~ 4 . The results in the last two lines are obtained 
for (Im A+/ Im A~)Bt/(Bg\/Rs) = ±1. For an explanation of the * see caption of 
Table §. 

6.3.3 Scenario C 

Within this scenario it is possible, in principle, to have a partial cancellation of the 
SUSY contributions to e' / e generated by Z-penguin and chromomagnetic operators. 
Given the strong RGE bound (|87|) , this possibility has only a minor impact on the 
upper bounds of both BR(A^ — ► tt°uu) and BK(K + — > tt + vv), with respect to 
scenario B. The only difference is that a sizable enhancement can also occur for 
ImA t > 0, if _R s ImA~ is positive and compensate for the negative contribution to 
{e'/e) generated by the Z penguin. This would then allow large values of K — > tiuu 
widths also within scenario I. This case is shown in Table As can be noticed, the 
upper bounds for the two neutrino modes within scenario II and III are the same as 
in Table ^ (with Re(e'/e) < 30.4 • 10~ 4 ) but, as anticipated, sizable enhancements 
occur also within scenario I. Due to the additional independent SUSY contribution 
to e'/e, in all cases (I-III) the upper bounds of K — > -nvv widths are insensitive to 
the experimental constraints on e'/e and depend only on the maximal value of At. 

More interesting is the case of Kl — ► ir°e + e~, sensitive to both ImAf and the 
SUSY contribution to magnetic operators. Also in this mode the largest enhance- 
ments occur when both ImAt arid R g lmA~ are positive, so that |i? s ImA~| can 
reach its maximum value. As shown in Table |6|, in this case one can reach val- 
ues of BK(Kl — > 7r°e + e - )dir larger than in scenarios A and B. An evidence of 
BR(A^ — > 7r°e + e~)dir ~ 10~ 10 would provide a clear signature of this particular 
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(though improbable) configuration. 

We finally note that, within scenario C, by relaxing the RGE bound (|7|) it is 
possible to recover the maximal enhancements for the rare decays pointed out in 
| |16| |. Needless to say, this possibility is rather remote, as it requires a few fine- 
tuning adjustments. However it is interesting to note that in the near future it 
could be excluded in a truly model-independent way by more stringent bounds on 
BR(i^ + — > 7i + uu). Indeed if BR,(Kl — > 7r°z/z/) > 2 • 10 -9 one expects from isospin 
analysis |33|] that BR(i^ + — > n + uu) > 4.6 • 1CT 10 , not far from the recent upper 



bound on this mode obtained by BNL-E787 34 



7 Summary 

In this paper we have analyzed the rare kaon decays Kl — > 7r°z/z/, K + — > 7i + uu, 
Kl —>■ ir e + e~ and the CP violating ratio e'/e in a general class of supersymmetric 
models. We have argued that only dimension-4 and 5 operators may escape the phe- 
nomenological bounds coming from AS = 2 transitions and contribute substantially 
to AS* = 1 amplitudes. On this basis we have introduced three effective couplings 
which characterize these supersymmetric contributions: A t for the Z penguin and 
for the magnetic ones. lmA t enters all rare decays and e'/e, ImA~ only e'/e while 
ImA+ only Kl 7c°e + e~. KeA t is important for K + — > 7i + uu and Kl — > fi + . 
Since ImA~ and ImA+ are expected to be similar in magnitude, a connection be- 
tween e'/e and Kl — > 7c°e + e~ follows in models with small ImA t . 
We have demonstrated explicitly that 

• the size of ImA^ is dominantly restricted by the present experimental range 
of e'/e; 

• the size of lmA t > is bounded by the minimal value of e'/e; 

• the size of Im At < is bounded by the renormalization group analysis (RGE) 
combined with the experimental values on ex and AM^; 

• the size of Re A t is bounded by Kl — > ^ \i~ and RGE. 

The imposition of the RGE bounds on the effective couplings has a considerable 
impact on the upper bounds of rare kaon decays (e.g. compare Table [| to Tables || 
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and so that the maximal branching ratios are found to be substantially lower than 



those obtained in [L5| . Given the important role of this bound it is worth empha- 
sizing that it requires more theoretical input than the low-energy phenomenological 
bounds usually taken into account within the mass-insertion approximation. Indeed 
it requires a control on the degrees of freedom of the theory up to scales of the order 
of 10 16 GeV. 

In order to accurately describe the relations between e' / e and the rare decays we 
have performed a numerical analysis in three basic scenarios: 

Scenario A: [ImA t = 0, ImA± ^ 0]. 

Scenario B: [ImAj ^ 0, ImA^ = 0]. 

Scenario C: [ImA t ^ 0, ImAj ^ 0]. 

In each of these scenarios we have considered three scenarios for the CKM factor At: 

Scenario I: X t is taken from the standard analysis of the unitarity triangle. 

Scenario II: Im X t = and Re X t is varied in the full range consistent with the 
unitarity of the CKM matrix. 

Scenario III: At is varied in the full range consistent with the unitarity of the 
CKM matrix. 

As we have discussed, scenario A with scenarios I or II for the CKM matrix is 
most natural within supersymmetric models with approximate flavour symmetries. 
However the other scenarios cannot be excluded at present and we have analyzed 
them in detail. Our main findings, collected in Tables [l|-|4] and |] are as follows: 

• In scenario A there is room for enhancement of BR(A^ ~~ > 7r°e + e~)dir by up to 
one order of magnitude and of BK(K + — > n + uu) by factors 2-3 over the Stan- 
dard Model expectations. BR(A^ — > remains generally in the ball park 
of the Standard Model expectations except for scenario II, where it becomes 
vanishingly small. 

• In scenario B, with the Standard Model values of ImAt (I), enhancements 
of BK(Kl — > ti^vv) by factors 2-3 and of BR(Kl — > 7r°e + e _ )di r by factors 
3-5 are still possible, while BR(A' + — > 7i + uu) can be enhanced by at most 
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a factor of 2. On the other hand, in scenarios II and III enhancements of 
BR(K L — > 7r°e + e~) dir and BR(K L — > ix^vu) by one order of magnitude and of 
BR(K + — > 7T + z/z/) up to a factor of 3 over Standard Model expectations are 
possible. These upper limits are dictated by the RGE bounds. 

• In scenario C enhancements of rare-decay branching ratios larger than in sce- 
narios A and B are only possible if ImA~ and Im A ( have the same sign so 
that the contributions of the chromomagnetic penguin and Z°-penguin to e' / e 
cancel each other to some extent. As a consequence the restrictions from 
e'/e are substantially weakened and what matters are the RGE constraints. 
In this rather improbable scenario one order of magnitude enhancements of 
BR(K L — > rc^uu) are possible even if the standard determination of A t is valid 
and BR{Kl — > 7r°e + e~)di r could reach the 10~ 10 level. On the other hand 
BR(K + — > -K + vu), being mainly sensitive to ReA t and ReA f , stays always 
below 3 • 10~ 10 as in scenarios A and B. 

We observe certain patterns in each scenario which will allow to distinguish be- 
tween them, and possibly rule them out once data on rare decays and improved data 
and theory for e'/e will be available. In particular in the near future with more 
stringent bounds on BR(K + — > 7i + uu) the most optimistic enhancements (like those 
occurring in scenarios C or Bill) could be considerably constrained. In the more dis- 
tant future, a clean picture will emerge from the measurements of BR(Kl — > 7r°z/z/) 
and BR{K L -> 7r°e + e -) dir . 
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